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Abstract 



^ ' We solve the coordinate space Faddeev equations in the continuum. We 

C**"* , employ hyperspherical coordinates and provide analytical expressions allow- 

ing easy computation of the effective potentials at distances much larger than 
the ranges of the interactions where only s-waves in the different Jacobi co- 
ordinates couple. Realistic computations are carried out for the Borromean 
halo nuclei 6 He (n+n+a) for J" = ± ,1 ± ,2 ± and n Li (n+n+ 9 Li) for |" 



2 > 2 > 

, | . Ground state properties, strength functions, Coulomb dissociation cross 

sections, phase shifts, complex S-matrix poles are computed and compared to 
available experimental data. We find enhancements of the strength functions 
at low energies and a number of low-lying S-matrix poles. 
PACS numbers: 21.45.+V, ll.80.Jy, 21.10.Dr, 21.60.Gx 

1 Introduction 

■ The present paper is part of a sequence discussing the general properties of three- 

body halos jl], ||, |[ H ||, ||]. These papers all deal with three-body systems, weakly 
bound and spatially extended compared to the energy and range of the two-body 
interactions j?], ||, ||, [k|. Most of the detailed information about halo nuclei is 
obtained from reaction experiments jy], [L^ |l3[ [l^] . Fragmentation reactions of 
three-body halo nuclei were studied in the sudden approximation for 11 Li and 6 He 
with special emphasis on the effects of final state interactions, which in other words 
are the effects of the two-body continuum ||, [l?], [ll| . 

Borromean systems, where no binary subsystem is bound, are particularly inter- 
esting three-body halo candidates. They have by definition a relatively low binding 
energy. The established nuclear prototypes are 6 He (n+n+a) and n Li (n+n+ 9 Li) 
and other examples are expected further up along the neutron dripline. The struc- 
ture of 6 He is fairly well understood whereas the structure of 11 Li is still controver- 
sial. The reason is essentially the large amount of knowledge, respectively the lack 
of knowledge, about the two-body subsystems. 

The number of bound states for Borromean systems is almost always limited 
to the ground state. The effective two-body interactions must be weak enough to 
exclude bound states and strong enough to bind the three-body system. Therefore 
one or more two-body resonances must be present at low energy. Then the low- 
lying continuum three-body spectrum would inevitably have rather complicated 
structure. This has strong implications for the analyses and the understanding of 
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the data accumulating from experiments with nuclear halos and Borromean nuclei 

The three-body continuum problem has been the subject of numerous investi- 
gations and tremendous progress has been achieved in recent years j27| ]. However, 
a number of problems still remain and unfortunately they include reliable com- 
putations of Borromean continuum spectra p8[. Other difficult problems include 
the continuum three-body Coulomb problem [[29| and scattering above threshold 
of one particle on a two-body bound structure as for example nucleon scattering 
on deuterons (3(], |3l| . The problems can crudely be divided into those dealing 
with the more technical aspects like the accuracies of the employed analytical and 
numerical methods, and physics issues addressing the behavior of the binary sub- 
systems, especially knowledge of the effective two-body interactions. It is necessary, 
but not always easy, to distinguish between these difficulties. Different treatments 
are usually needed for short-range and long-range interactions and for energies be- 
low or above possible two-body thresholds. The concepts and problems described 
above are to a large extent general and of interest in other subfields of physics 
H,|||||34H||6H3. 



The general discussions of structure and break-up reactions of halo nuclei should 
be extended to the three-body continuum. Specific: investigations are available, but 
imprecise |2C], |l| H) ■ ^he technical difficulties formulated in coordinate space 
are related to the necessary computation of the behavior of the effective potential 
at large-distance. Fortunately, a method treating the large distances analytically 
and the short distances numerically has recently become available (3^, pSj . The 
method is very powerful in structure calculations as demonstrated by the succesful 
investigation of the Efimov effect ]8|, |3^, |4C , 41 1 . However, the implementation has so 
far essentially concentrated on bound structures, but generalization to applications 
in the three-body continuum is straightforward. 

The purpose of this paper is to (i) describe the details of a method to compute 
low-energy three-body continuum spectra for particles with or without intrinsic 
spin, (ii) derive asymptotic large-distance expressions allowing simple computations 
of the corresponding effective three-body potential for arbitrary angular momenta 
and arbitrary short-range two-body potentials, (iii) apply the method in detailed 
realistic numerical computations of the continuum structure and various observables 
for 6 He and n Li. 

The paper generalizes first the analytic results obtained for s-waves and square 
well potentials jy] . Then the method is applied to detailed studies of the continuum 
structure for the Borromean nuclei 6 He (n+n+ 4 He) and u Li (n- 
reports describing some numerical results are available in the literature |26 

After the introduction we give in section 2 a general description of the method 
and then we concentrate on two cases of special interest. In section 3 and 4 we 
compute in details the properties of 6 He and 11 Li, respectively. In section 5 we 
give a brief summary and the conclusions. A convenient general expression for 
the transformation between different sets of Jacoby coordinates is derived in the 
appendix. 



2 Theory 

We shall consider a system of three interacting inert "particles". Their intrinsic 
degrees of freedom are frozen and only the three-body (halo) degrees of freedom 
shall be treated here. In this section we first describe the general method of solving 
the Faddeev equations using hyperspherical coordinates. In particular we specify the 
boundary conditions at large distances by introducing the S-matrix or equivalently 
the i?-matrix. The angular equations at large distances are then treated essentially 
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analytically. We then consider a system of two identical neutrons surrounding a 
core first with finite spin and then with spin zero. Finally in this section we give 
the expressions for strength functions and Coulomb dissociation cross sections. We 
shall follow the method and the notation established previously (?], ||, |j| |8[ gjj, |2(| . 



2.1 Method 

The k'th particle has mass to^, charge eZ^, coordinate rj, and spin s^. The two- 
body interactions between the particles i and j are . We shall use the three sets 
of Jacobi coordinates (x^y^) and the corresponding three sets of hyperspherical 
coordinates (p, Qi)=(p, ai, ^yj, see for example || |7j. The volume element 
in terms of one of the sets of hyperspherical coordinates is given by p 5 dfldp, where 
dr2 = sin 2 a cos 2 adadQ x d£l y . The kinetic energy operator is 

/ £)2 K A 2 \ 

(1) 

sin(2a) da 2 v ' ' sin 2 a ' cos 2 a ' 

where the angular momentum operators I 2 and I 2 are related to the x and y degrees 
of freedom and m is a normalization mass arising from the definition of p. In the 
following m is assumed to be the nucleon mass. 

The total wave function \P jm of the three-body system (with total spin J and 
projection M) is written as a sum of three components V'jm > which in turn for each 
p are expanded on a complete set of generalized angular functions &n\p,tti) 




3 oc 3 

*jm = J2^M(xi,yd = ^72 E^frOE^fonO , (3) 

i—l ^ n—1 i—1 



where p~ 5 ^ 2 is the radial phase space factor. 

The angular functions are now for each p chosen as the eigenfunctions of the 
angular part of the Faddeev equations: 

~ (A 2 - A n (p)) *W + V$ fc (*W + + *<*>) = , (4) 

where {i,j, k} is a cyclic permutation of {1, 2, 3}. 

The radial expansion coefficients f n (p) are obtained from a coupled set of "ra- 
dial" differential equations 0, |42|] , i.e. 

d 2 2m(E-V 3 (p)) , 1 A 15\ \ 

A„(/9) + — - Qnn MP) 



9p 2 a 2 V V 4 

d_ 
"dp 



hPun'-^ + Qnu') fn'(p) , (5) 



where E is the three-body energy, Vs(p) is an anticipated additional three-body 
potential and the functions P and Q are defined as angular integrals: 

p nn ,( P ) = E / dn^fafijJU^fon) , (6) 
Q„ n ,( P ) = E / dn*J? , (p,fi)^^ ) (p J n) ■ (7) 
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The diagonal part of the effective potential is then 



ft 2 , , 15, _ , 



. )/;; . (A„(P) + -j)p- 2 ~ Qnnj + V 3 ( P ) . (8) 

For Borromcan systems the coupling terms P and Q approach zero at least as 
fast as p -3 j?], |38|, We can then choose those solutions to eq.([|) where the 
large-distance (p — » oo) boundary conditions for /„ ' are given by 

fk n ' ] (p) - *n,n'^ -) («P) ~ S n , n ,F^ (K P ) . (9) 

The S- matrix introduced here is a unitary matrix, n 2 = 2mE/h 2 and F^' are 
related to the Hankel functions of integer order by 



±Z K p±y(X„ + ^) 



(10) 



where K n is the hyperspherical quantum number corresponding to the value K n x 
(K n + 4) approached at large distance by the angular eigenvalue A„ . The continuum 
wave functions ^ jm are orthogonal and normalized to delta functions in energy. 
Sometimes it is more convenient to work with the f?-matrix given as R = i(l — 
S)/(l + S). The boundary conditions must then be changed into sin and cos instead 
of the exponentials in eq. (|I(]) . 

By diagonalization of the S- (or i?)-matrix we obtain eigenfunctions and eigen- 
phases. These phase shifts reveal the continuum structure of the system. In par- 
ticular, a rapid variation with energy indicates a resonance. A precise computa- 
tion of resonances and related widths can be done by use of the complex energy 
method, where eq.(^) is solved for E — E r — iT/2 with the boundary condition 

fn oc J Hr;} +2 ( k p) ■ These solutions correspond to poles of the S- matrix p3| . 
2.2 Angular eigenvalue equation 

The angular functions <I>^(p, fli) are expanded in products of the three-body spin 
functions xi*\ Sm anc ^ spherical harmonics Y£^ mx (£l Xi ) and Yg yrriy (Q yi ). The or- 
bital angular momenta and their projections associated with x and y are (£ xi m x ) 
and (£ yi m y ) while the spins of the two particles connected by the x coordinate cou- 
ple to the spin s x , which coupled to the spin s y of the third particle results in the 
total spin S and its projection m s . Indicating these angular momentum couplings 
by <g) the result can be written 

= E ^a;^ 1 [y*7^M® x z ySm ,] JM , (id 

where sin(2o!i) is a factor related to phase space and 

Y t L f u L (Q Xi ,il yi ) = [Y^ mx (n Xi )®Ye ymy (n yt )] LML , (12) 

where the projections of the intermediate couplings are given although the final 
result is independent of them. 

To solve the angular Faddeev equations the components in eq.Q must be 
expressed in one Jacobi coordinate set, say labeled by i. The wave functions 
4>ni i Ls s (p, OLj)j sin(2aj), which only depend on cij and p, are first expressed 
in terms of the i'th set of hyperspherical coordinates. The equations are multiplied 
from the left by the square bracket in eq.(fTl) and subsequently integrated over 
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the four angular variables describing the directions of Xj and yj. The interaction 
Vjk, only depending on the distance between the particles, is independent of these 
angles. 

The operator describing this transformation from the j'th to the i'th Jacobi 
coordinate system is denoted Rij. This operation maintains both total spin and 
total orbital angular momentum. The result of the transformation from a specific 
set of angular momentum states i' x i' y L is projected on the set £ x £ y L. This operator 

t l l' l' L 
R y * y is then given by 



sin(2a,) 



dfL.dfL 



Y, 



lm l 



iS^Xi 1 ^"iji ) 



When the two-body interaction is assumed to be diagonal in the total two-body 
spin we now rewrite the angular eigenvalue equation in eq. (Q) as 

\ d£f + sin 2 ai + cos 2 a 4 + P V i (P sma *) ~ v n\P) ) Pnl x £ y L S!e SKP' a ^> 



+p*sm(2a t )v:* s (psma l )j: e , e , s , [ C^^f* 



sin(2ctj ) 



^^s^s^'S^ik 



, , . (p.OCi.) 

! ~ e y Ls L S 

n(2a k ) 



o, 



(14) 



where i%(p) = A„(p) + 4 and the reduced and spin averaged interactions are given 
by ^ 

«i" S (aO = {Xs xSv Sm s \^2 V ik{— )\Xsls v Sm s ) ( 15 ) 
h pjk 



with mHj k = mjirLk/(mj + rrik). The coefficients C™ ax , s , expressing the overlap of 
the spin functions, are defined by 



s x s v Sm s I As' s^,Sm 



(k) 



(16) 



These matrix elements are independent of m s , symmetric, i.e. C^", s = C^f sS and 
diagonal in s x and s' x for i — k, i.e. C^, s = <J SS '. 



2.3 Large-distance angular eigenvalues 



For large p only small a-values contribute in eq.(14) to the terms proportional 
to p 2 v1* s (psincti). These potentials are assumed to have short ranges and they 
vanish consequently for large p for all ai except in a narrow region around zero. 
We assume that they vanish exponentially or at least as fast as 1/p 3 , for distances 
larger than the ranges of the potentials 0, |dj . The two terms described by the 
transformations Rij and Rik in eq. (^) can then be approximated by their expansion 
to leading order in the variable ai. 

We show in appendix A that all partial waves decouple to leading order in a, 
or in 1/p with the essential exception of s-waves in the rr-degree of freedom, i.e. 
the components with i Xi = 0,i — 1,2,3 and the total orbital angular momentum 
L — £ yi . This means that an expansion in powers of ai of the terms obtained from 



■5 



the transformation R 
for 4 = 4 = 0, € 



y only provides non-zero contributions in the limit of ctj = 
L. These finite contributions are for a,; = found to be 



v 



b: 



OLOLL 



<>nOLLs> x s(P> a 3) 

sin (2a j) 



= (-1) 



sin(2(^ fe ) 



tant^fc = 



lmk(rrii 



- m k 



mirrij 



(17) 



(18) 



where {i,j, k} again must be a permutation of {1,2,3}. Higher order terms in a, 
are neglected. Non-zero 4-values had produced leading terms of at least first order 
in ai in the expression analogous to eq.(fl7|). 

Thus for non-zero t ^-values the angular eigenvalue equations in eq. ( |li| ) decouple 
asymptotically and reduce to 

d 2 
da 2 



+p 2 v^ (psinaj 



sm ai 



Ls x S 



{P,Oii) 







(19) 



for all sets of values of 

,2 



7^ 0, £ y , L, s x , S. The large-distance asymptotic eigenvalues 
v~ = (K + 2) 2 related to these partial waves approach the hyperspherical spectrum, 
where K is odd or even natural numbers depending on the parity. This asymptotic 
behavior is reached on a distance scale defined by the short range of the interactions 
v°* S in eq.©. 

For 4 = insertion of eq.(|l7]) into eq. (|l4|) gives instead the three coupled 
asymptotic angular equations 

da 2 



2 +l${p,OCi) ) ^oii S „s(P)«i) 



2a. l ^l) L C { [> SxS p 2 vr S {psina l ) 



(20) 



C 



(i) 

Ls x S 



K 2 (p, a t 



E & L 



L(L+1) 



+ p 2 v^ s (pain a t ) 



"UP) 



Polls' s(p. V*) 



sin(2</?* 



sm(2(fij) 



(21) 



(22) 



Also these eigenvalue solutions j/ 2 converge towards the hyperspherical spectrum as 
p increases. However, due to the coupling the asymptotic values are now approached 
over a distance defined by the scattering lengths, which might be very much larger 
than the ranges of the interactions. 

As mentioned above the potentials p 2 vf* s (psmai) vanish for large p for all 
ai except in a narrow region around zero. The conditions for the effective range 
approximation therefore become better and better fulfilled as p increases and any 
potential with the same scattering length and effective range would lead to the 
same results. Let us then in the region of large p use square well potentials 



Vjk(r) = — S , Q J ' l (sjSfe)6(r < i?* 3 " 5 ), or equivalently expressed by the reduced quan- 
tities v^ S (x) = -s^\s x S)Q{x < X-" S = R^ S fi jk ), where the range and depth 
parameters are adjusted to reproduce the two-body scattering lengths and effec- 
tive ranges of the initial potential. The corresponding solutions are then accurate 
approximations to our original problem at distances larger than 2R 



s x S 
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The square well potentials v? x (psino^) are zero in region II defined by a, > 



a[) l) (s x S) 



encsm(X^ s /p). Then eq.| 

d 2 L(L + 1) 
da 2 cos 2 on 



is especially simple, i.e. 



(23) 
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and the solutions, vanishing at a^ = 7r/2, are given by 



<t>nOLL Sx s(P> a ) = A nOLL Sx S P L{Vn,a) 



Pi(i/„, a) = cos L a 



for arbitrary constants -A^QLLs s 



d_ 1_ 

da cos i 



("-?) 



(24) 
(25) 



The potentials v^" S (psinai) are finite and constant for large p in region I defined 
by ai < Oq (s x S) <C 1. Then eq.(20) is approximately 



+ = 0)) = -2a i (-l) i p 2 S «( Sx -5)C 



(i) 



(26) 



where the wave functions in C^ s s in eq.(|22]) must be (f'nOLLs s- The solutions to 
eq. (|26|) , vanishing for aj = 0, are then 

& LSx s(P: «) = ^ll,,s sin(a Kl (p, a = 0)) - 2a(-l)^ ^^^ C^ s (27) 
for arbitrary constants B^q LLs s , where /tj are defined in cq. (pf|) . 



Matching the solutions, eqs.(|24|) and (|27|), and their derivatives at a% — a^ (s x S) 
gives a linear set of equations for A]lq LLs s and bQ ll s with given L and S for 
i = f,2,3 and all possible s x . Physical solutions are then only obtained when the 
corresponding determinant is zero. This is the quantization condition for v 2 (or 
A) and as such the eigenvalue equation determining the large-distance asymptotic 
behavior of X(p). 

The square well solution in eq.(|27]) is not exact since the first order expansion 
in ai is used in eq.(pOf) and (^6|) and consequently also in the last term of eq. (|27|) . 
Improvements could be obtained by using eq.(24) in eq.(A|l5|) and thereby changing 
the right hand sides of eqs.(|l7|), (|20|), ( |2^ ) and (p7|). For L = these expressions 
are given in 

Also the eigenvalue equation for non-zero ^-values in eq.([l9]) can be solved for 
square well potentials. For ai > Oq (s x S) in region II, we have the equation 



da 



4(4 + 1) 4(4 + 1) 2 \ (i) 

1 1 1.2 „.. + c J 2a . " "n J <fHt y L. m s(P> «0 = (28) 



sin ai 

and the solutions vanishing at 7r/2 are then 



N„it sin** a cos** aP^ +1 / 2 ^ +1 / 2 (cos(2a)) (29) 



i(0 



where p^+VVy+i/s are th e j aco bi polynomials 



for arbitrary constants ^ £ S:r s, 
and N n i x i y are normalization constants given explicitly in appendix A. 

For ai < aQ\s x S) <C 1 in region I eq.(|i9|) is approximately 



d 2 | 4(4 + 1) 

da 2 a 2 



with the solutions vanishing as a* at a, = 



^Li^s^") = B nle v L S!C sie*( aK i(P> a = °)) 



«J(p, a = o) = - 4(4 + 1) - P 2 4 l) («x5) - ^ 



(30) 

(31) 
(32) 
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where B^j e Ls s is an arbitrary constant and jt x is the spherical Bessel function, 
i.e. the usual solution to the radial two-body Schrodinger equation for an angular 
momentum £ = £ x . 

Matching the logarithmic derivatives of the solutions in eqs.(p9|) and ( j3l|) then 
provides the quantization condition for v 2 (or A) and therefore the large-distance 
asymptotic behavior of A(p). 



2.4 Large-distance behavior for two neutrons and a core 

We shall now consider a system of two neutrons (labeled 2,3) and a core (labeled 1) 
with spin s c . For a given total spin and for each set of orbital quantum numbers the 
six possible components, ^ng^e Ls^S' * = 1> ^ 3, each with two s^-values, are related 

to the three-body spin wave functions xi*=o ij X^Ls ±1/2 > X, =« ±1/2' wnere the 
first set of Jacobi coordinates corresponds to the ^-coordinate between the two 
neutrons. Due to the Pauli principle only three of these wave functions (f> are 
independent and the remaining components are determined by antisymmetry, i.e. 

4>nl xf > yLs *S = for odd 4 + s x , (33) 

Specifically the three independent s-wave components can be characterized by 
L, S and one value of s x , i.e. ^ = tf^ LLogt <j)f s = (/> ( ^ LLsc _ 1/2S , 4>f s = 
(3) 

0nOLLs +1/2S- These components are coupled over a distance defined by the scat- 
tering lengths whereas all other partial waves decouple for large p above a distance 
scale defined by the range of the interactions. 

(i) 

The components q> LS obey for large p the coupled angular Faddeev equations in 
eq. (|20|) , where the coefficients cj^ s = C^l s (s Xl = 0, s X2 — s c — i, s X3 = s c + |) 
explicitly are given by 

^(1) _ 9^12 ^j 2 s(y) _ 9^12 <f>LS W (oz\ 

°LS — ZO 0,s c -l/2,S sin(2 V ) ZL/ 0,s c + l/2,S sin(2 V ) y 00 ) 

ri 2 ) — r 12 <t>Ls(v) , r 23 tSUfl . r 23 ^i 3 s(y) 

^LS ~ K - y 0,s c -l/2,S sin(2 V ) ^ s c -l/2,s c -l/2,S sin(2£) L/ s c -l/2,s c + l/2,S sin(2£) 

<o(3) _ ^13 ^a(y) , ,-23 j£|(g) _ /-i23 ^ls^) 

^LS Ly 0,s c + 1/2,S sm(2 v ) ^ Ly s c -l/2,s c + l/2,S sin(2£) °s +l/2,s c +l/2,S sin(2y) ' 

where we omitted the argument p in the functions and further defined ip = if2 — 
tp>3, <p = (p±. For S = s c all terms are present, but for S = s c ± 1 the first term in 

(i) 

C LS should be removed together with the equation corresponding to i = 1. 
The spin overlap coefficients are explicitly given by 



Q) 2 s c -l/2,s c — Co 3 < ic -l/2, ;ic — J 2 S _|_ I 



Co,s c + l/2,s c — Qm c + 1/2,s c — Y 2 S +1 



c: 

^23 

u s c -l/2,s c -l/2,s c +l 



2.; 





^23 V'M^ + 
°«.-l/2.-.+l/2,«. - 2Sc + l 


1) 


(38) 


i 

-1/2, 


^23 ' 

So -l/2,s c — Ly s c + l/2,s c + l/2,s c — 2 Sc -| 


- f 


(39) 


,-23 


✓~<23 

-l/2,s c + l/2,s c -l — Ly s c -l/2,s c + l/2,s c ±l " 


= 


(40) 


c 2 ^ 


Ai23 

l/2,s c -l/2,s c -l — °s c + l/2,s c + l/2,s c + l " 


= 1 


(41) 
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The potentials p 2 Vi(p sin a) approach for sufficiently large p the zero-range po- 
tentials, where the sensitivity to the shape disappears. Any potential with the same 
scattering length and effective range would then lead to results accurate to the order 
p~ 2 . We shall therefore for convenience use such equivalent square well potentials, 
where the solutions to eq.(pp|) then again are given by eqs.(p4|) and (pTj), and the 
large-distance physical solutions are obtained as described above. 

We shall now consider a system of two neutrons and a core with spin s c = 0. 
All quantities where s c — \ appears as an index should now be substituted by zero. 
Then the three coupled s-wave equations in eq.(po|) reduce for S — to two as seen 
from cq.(|3^) where C^g then is zero. These equations are to leading order in a 
(large p) explicitly given by 

9 2 L{L + 1) 2 . . , a \ 
^2 + ~S " + Ann PSIeqj - v 2 42 

)OL\ COS^ Ct\ J 

x^Vai) = -2a 1 (-l) L j o 2 « NN (psina 1 )C[ 1) , 

p 2 v Nc (psma 3 ) - v 2 I (43) 



d 2 | L(L + 1 ) , . , ; , 

dot\ COS 2 a 3 

x^ 3) (p,« 3 ) = -2a 3 (-l) i j0 2 UNc (psina 3 )ci 3) 
where v NN (xt) = v[ 00 \xx), v- Nc (x 2 ) = v { 2 1/2,Q) (x 2 ) = v^ /2fi \x 2 ), = </>£o> c 



C« and 



W 
L 



L sm{2ip) ' sin(2<£>) sm(2£) ' 



The equivalent square well solutions are again given by eqs.(|2J) and (^7j), and the 
large-distance asymptotic behavior are obtained as described above. 

For S — 1 only eq. ([43]) remains for s-waves now with — 4>^ (if)/ sin(2y>) . 
The square well solution and the large-distance behavior are then easily obtained. 

2.5 Strength functions and Coulomb cross sections 

The strength functions d ^J A describing electric multipole excitations of the ground 
state \Jq°) into the continuum state \nJ n E) are defined by 

dB EX (E) = 

9,.Tn -4- 1 

nj 71 



dE 2J, 



-^^KnJ^HM^A)!!^)! 2 , (45) 



M(EX, fi) = J2 eZ k r%Y Xll (r k ) ■ (46) 

k=l 

in terms of the reduced matrix element and the electrical multipole operator M. 
The corresponding sum rule is 

/ di ^l^ = ^E^ 2 W°K 2A |^> , (47) 

where only the core contributes for a system of two neutrons around a core. 

Nuclear excitations of monopole type are possible with the corresponding oper- 
ator p 2 = J2k( r k ~ R-c) 2 , where R c is the coordinate of the center of mass. The 
related strength function d ^J and the sum rule are then 



dE 

dN E0 (E) 

dE 2J 



L-J2\(nJS"E\p 2 \J^)\ 2 , (48) 



d£^f^ = ^ (W°IpVo°> ~ (Wm 2 ) , (49) 
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The Coulomb dissociation cross section can now be computed in the high beam 
energy limit where the approximation of straight-line trajectories is valid and only 
one photon is exchanged between projectile and target. The cross section is then 
obtained by multiplying the electromagnetic transition matrix elements dB ^ E ^ 
from ground state to continuum states with the virtual photon spectrum nE\(w), 
which is given by G3 



n_B 2 (w) 



n E \{uj) 
2 



2 r 2 

Z rj2 L 



2 C 4 

- z t a ~z 



2 1 



2(2 
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and the resulting differential cross section is 
da E i(E) _ 167r 3 a 



dE 

daE2(E) An 3 a 



dE 
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9 

E* 
ft c 



n E1 (E /%) — 



n E 2(E*/h) 



e* dE* 
1 dB E2 {E* 



dE* 



(51) 

(52) 
(53) 



where K is the modified Bessel function, a = e 2 /he, £ = v is the beam velocity, 

7 = l/yl v 2 /c 2 , Z t is the charge of the target and E* = huj = Ef — Ei, where the 
final and the initial energies are labeled by f and i. The dipole is usually by far the 
largest contributor. In any case for most halo nucei (the quadrupole excitation for 
6 He is an exception) the information about higher-lying angular momentum states 
is not experimentally available and very difficult to predict theoretically due to the 
lack of knowledge about the binary subsystems. 



3 The 6 He-system as n+n+ai 

he three-body model developed above can be tested on the 6 He-system, which has 
been studied as the simplest prototype of a halo nucleus. The advantage is that the 
details of the low-energy two-body interactions are very well known experimentally 
and the particles only have high-lying excited states. The resulting three-body 
properties are therefore much less uncertain and related to the technique rather 
than to the lack of information about the subsystems. In this section we shall first 
study the influence of the remaining uncertainties in the model, then predict physical 
properties of the three-body system and along the way compare with available data. 



3.1 Interactions and numerical details 

We consider 6 He as two neutrons and an inert 4 He-core. The two-body interactions 
should in principle only reproduce the low-energy scattering data which exclusively 
influence the computations of spatially extended halo systems. Except for very 
accurately needed details it is even quantitatively sufficient to reproduce the scat- 
tering lengths and the effective ranges of the appropriate partial waves. This initial 
conjecture p5| is now confirmed and details of the short-range behavior of the two- 
body interactions are not needed ^ ^, [l0|, |4lJ] . We shall therefore essentially always 
maintain the same radial shapes of the interactions. 

The neutron-neutron interaction reproduces the low-energy properties of free 
nucleon-nucleon scattering. We have tried several parametrizations, i.e. the sim- 



ple neutron-neutron s-wave potential, -31 MeV exp(— r 2 /(1.8fm) 2 ), from |45|, the 
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Figure 1: The computed s-, p- and rf-wave neutron-a phase shifts (solid curves) 
compared to the values (triangles and circles) extracted from scattering experiments 



1471 148| . The interactions are given in eq.(54) 



extension to other partial waves in |J, the accurately adjusted nucleon-nucleon po- 
tential from |l(J and previously known standard potentials as that of f4£)| . The 
three-body results can hardly be distinguished from each other and we shall here 
only present results with the interaction from [0]. 

The neutron-a interaction is parametrized to reproduce accurately the s-, p- and 
d-phase shifts up to 20 MeV. We use again gaussians for the radial shape and allow 
an ^-dependence of strengths and ranges, i.e. 

pKf = °) = 48.00 exp(-r 2 /2.33 2 ) 
V 7 i e c =1) = -47.40 exp(-r 2 /2.30 2 ) - 25.49£ • s„ exp(-r 2 /1.72 2 ) (54) 
V^= 2) = -21.93 exp(-r 2 /2.03 2 ) - 25.49£ • s„ exp(-r 2 /1.72 2 ) , 

where the strengths are in MeV, the lengths are in fm, s n is the neutron spin 
and I is the relative orbital angular momentum. The repulsive s-wave potential 
corresponds to a scattering length of -2.13 fm and an effective range of 1.38 fm. 
The energies and widths of the p-resonances defined as poles of the S'-matrix are 
E(p 3/2 ) = 0.77 MeV, T(p 3/2 ) = 0.64 MeV, E{p 1/2 ) = 1.97 MeV and T{p 1/2 ) = 5.22 
MeV, respectively. The phase shifts from this potential are in Fig. 1 compared with 



the results obtained from scattering experiments |£7j, |48 |. 

Other parametrizations are possible even for the same radial shape of the two- 
body potential. They differ in the number of two-body bound states of which the 
lowest s-state is occupied for 6 He by the core neutrons and therefore subsequently 
has to be excluded in the computation. For one bound s-state the interaction for 
£ = is 

F n (f=°) = -75.06 exp(-r 2 /1.53 2 ) (55) 

while the £=1,2 partial waves remain the same as in eq. (p4) . The s-wave scattering 
length and the effective range are the same as for eq.(F34|) although the potential 
now is attractive. 
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The three-body system computed from these two-body interactions is under- 
bound by about 500 keV. The required fine tuning is now obtained by adding a 
diagonal three-body force V3 (p) in eq. (||) . The idea of using the three-body force is 
to include effects beyond those accounted for by the two-body interactions. Thus 
two-body polarization effects are already included via the effective two-body inter- 
action. The remaining part must then involve all three particles simultaneously 
polarizing each other and therefore only effective at small p-values. We therefore 
use a three-body interaction only depending on p and still allowing for a depen- 
dence of the total angular momentum of the system. We tried both gaussian and 
exponential shapes, i.e. Vs(p) = S 3g exp(— p 2 /b 3g ) and Vs(p) — S 3e exp(— p/b 3e ). 

The range of the three-body force is by its definition related to the hyperradius. 
For 6 He, p=2 fm and 3 fm correspond roughly to configurations where the neutrons 
respectively are at the surface of the a— particle and outside the surface by an 
amount equal to their own radius. This distance can now be used directly as the 
range parameter or defined as the distance where the three-body potential assumes 
half of its central value. This means that b 3e = b p /\n2 and b 3g = 6 p /Vhi 2, where 
bp = 2 fm or 3 fm for the two different geometric configurations. 

The strength of the three-body interaction is finally for + adjusted to give the 
measured two-neutron separation energy 0.97 ± 0.04 MeV of 6 He. The different 
ranges and radial shapes has an influence on the spatial extension of the three-body 
system. For gaussian shapes we obtain root mean square radii of 2.45 fm for both 
the attractive s-wave potential with one bound state and the repulsive potential 
without bound states. The corresponding three-body interaction parameters are 
respectively b 3g = 2.9 fm, S 3g = -7.55 MeV and b 3g = 3.0 fm, S 3g = -3.35 MeV. 
For exponential shapes and repulsive s-wave potential, we obtain instead root mean 
square radii varying almost linearly from 2.61 fm to 2.56 fm for S 3e — —3.11 MeV, 
b 3e = 4.3 fm to S 3e = —4.77 MeV, b 3e — 3.0 fm. For 2 + we could instead fine tune 
to the well known resonance of energy 0.820 ± 0.025 MeV and width 0.113 ± 0.020 
MeV j49|. The three-body interaction parameters would then be J-dependent. For 
the repulsive potential we obtain b 3g = 2.061 fm and S 3g = —31 MeV for gaussian 
shapes. 

In the computations we include all possible s-, p- and c£-waves. We use a hyper- 
spherical basis for each of the Faddeev components with iV values up to about 150. 
The radial equations are integrated from zero up to p- values about 180 fm. Further 
arguments for these numerical choices can be found in |26]| . 



3.2 Solutions and iS-matrix poles 

The angular eigenvalues A„ are computed from eq. (|l4j) for total angular momentum 
and parity J 71 " = ± , 1 ± ,2 ± . These eigenvalues are closely related to the effective 
potentials in the radial equation eq. (j^) . Their large-distance behavior is essential 
and sometimes decisive as seen in the extreme case of Efimov states which owe 
their existence to a sufficiently negative value of A at very large distance [^8[ |ib]| . 
We show in Fig. 2 these spectra for the lowest spins and parities both computed 
numerically with the appropriate basis size and from the analytic expressions for 
the coupled s- waves. 

The asymptotic behavior obtained at large distance is in this case reached around 
40 fm. This does not mean that all interactions produce the same results at such 
a distance. It means that details are unimportant, but the scattering lengths are 
still crucial. A larger basis would have reproduced the results of the analytical 
calculations up to higher p-valucs. However, this is not needed, because we use 
the asymptotic solutions as soon as they are accurate enough. This improves both 
accuracy and speed of the computations. The finite size of the basis gives a too fast 
convergence to the hyperspherical spectrum. Without an independent calculation 
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Figure 2: The lowest angular eigenvalues A„ for 6 He (n+n+a) as functions of p 
for angular momentum and parity J* = 0^, 1 ± ,2 ± . The solid lines are computed 
by numerical integration and the dot-dashed lines are the large-distance asymptotic 
behavior obtained from eqs.([l2"|) and j43|). The neutron- neutron interaction is from 
|^| and the neutron-a interaction is from eq. (|54|) . Maximum if- values up to about 
150 are used in the basis. 



13 



40 




repulsive potential 
attractive potential 



30 



20 



Q. 




-10 







5 



10 



15 



P(fm) 



Figure 3: The lowest angular eigenvalues A„ for 6 He (n+n+a) as functions of p 
for angular momentum and parity J 71 — + for the repulsive (solid) and attractive 
(dashed) s-wave potentials in eqs.(|54|) and j55|). The neutron- neutron interaction 
is the same as in Fig. 2. The lowest diverging level for the attractive potential is 
omitted corresponding to our prescription for excluding the Pauli forbidden state. 
Maximum i^-values up to about 150 are used in the basis. 

it can therefore be difficult to assess the accuracy. 

The lowest level for each J 71 " usually contribute with the largest components of 
the wave function of both the possible bound state and the low-lying continuum 
states. For all six cases in Fig. 2 we find pockets which, except for the + ground 
state, are unable to bind the system, but still responsible for several low-lying S- 
matrix poles as we shall see later. 

The Pauli principle prohibits the core neutrons and the halo neutrons from occu- 
pying the same orbits. This has to be incorporated explicitly, since the three-body 
model only deals with particles without intrinsic degrees of freedom, except for their 
intrinsic spins. For a repulsive s-wave potential as described above, no bound state 
is present and no overlap has to be excluded. We also investigate another approx- 
imation where we use an attractive neutron-core potential with one bound s-state 
and the same scattering length and effective range. The lowest angular eigenvalue 
must then bend over and diverge parabolically towards — oo. This corresponds at 
large distances to a configuration where a neutron is bound in the doubly degenerate 
lowest s-state, which is Pauli forbidden for the halo neutrons. At smaller distances 
the probability, or the wave function, must be small, because otherwise a significant 
part of the halo wave function would be inside the core, the halo and core degrees of 
freedom would not separate and the three-body model would not be a good approx- 
imation. The effective potential at these small distances is then rather unimportant 
if the model is valid. Therefore a good and inexpensive approximation to include 
the Pauli principle is simply to omit the lowest diverging angular eigenvalue A from 
the computations ||. 

In Fig. 3 is shown the angular eigenvalue spectrum for J n = + for both the 
repulsive and the attractive s-wave potentials in eqs.(p3) and (|55|). The lowest 
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Table 1: The real and imaginary values (E r , V) (in MeV) of the two lowest S'-matrix 
poles E — E r — ir/2 for 6 He for various spins and parities. The interactions used 
in the upper part of the table are the same as in Fig. 4. The three-body interaction 
parameters are S 3g = -7.55 MeV, b 3g = 2.9 fm, S 3g = -31 MeV, b 3g = 2.061 fm, 
respectively for the first two and the last two columns. The 1~ poles in the middle 
are obtained with the same interactions except for an exponential shape for the 
three-body potential with parameters S 3e = —3.11 MeV, b 3e — 4.3 fm. In the lower 
part of the table the repulsive potential is substituted by the attractive potential in 
Eq.(|5^). The interactions used are otherwise unchanged, except for the three-body 
interaction parameters, where S 3g — —3.35 MeV, b 3g — 3.0 fm. The excitation 
energies are E* = E r + 0.95 MeV, E* = E r + 1.54 MeV for the left and right hand 
side of the table. 
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0.44 
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diverging level for the attractive potential originating from zero is removed from 
the figure as well as from the subsequent computations. The second level for the 
attractive potential, originating from 12, is almost identical to the lowest level from 
the repulsive potential from about p = 1 fm. The levels from these two potentials 
are remarkably similar even at smaller distances and they are completely identical 
in the large-distance asymptotic region. (Note that the figure only shows results up 
to 15 fm, where differences still can be seen.) One level originating from 32 must 
cross an empty region, and therefore deviate somewhat from all other levels, until it 
catches up with one of the levels from the other potential approaching 12 for large 
p. However, the lowest A-value(s) is dominating in the wave functions of interest 
here and we therefore should focus on the corresponding effective potentials. The 
differences in these potentials are small, but in precise computations they must be 
compensated in one way or another. Fortunately the means for such fine tuning is 
already present as a three-body potential, which is different for different two-body 
potentials. 

For the lowest spins and parities we give in Table 1 the lowest resonance energies 
and the related widths obtained as ^-matrix poles by the complex energy method. 
The different three-body forces in Table 1 can be considered to give the realistic 
range of the possible variation in the present model. In contrast, all previous com- 
putations did not produce three-body resonances in this low-energy region, except 
the established 2 + -resonance [^3[ |24|, |25| , For 0~, 1 + and 2~ we obtain identical 
poles, since the two three-body interactions only contribute at small distances where 
the effective two-body potentials completely dominate. On the other hand, for + , 
1~ and 2 + we find differences of up to 0.3 MeV and 0.18 MeV for the position and 
the width, respectively. The systematic shifts of the positions in the right hand side 
of the table arise due to the slightly different 2 + energy obtained by adjusting the 
parameters. 

The radial shape of the three-body force and the change from the repulsive to 
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Figure 4: The total effective diagonal radial potentials for 6 He (n+n+a) defined in 
eq.(||) (solid curves) as functions of p corresponding to the three lowest A's for ± , ^ 
and 2 ± . The dashed curves are the part remaining after removal of the generalized 
centrifugal terms, i.e. fi 2 {K + 3/2) (K + 5/2)/(2mp 2 ), where K(K + 4) = A(p = 
oo ) is the corresponding asymptotic hyperspherical eigenvalue. The interactions 
are the same as in Fig. 2 with an additional diagonal three-body interaction, i.e. 
Sz g exp(— p 2 /b\ g ) with Sz g = —31 MeV, bz g — 2.061 fm, added in all partial waves, 
except for J = + where S^ g = —7.55 MeV, bs g = 2.9 fm. The insets show the 
details of the lowest potentials. 
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the attractive s-wave potential both have an effect on the lowest S-matrix poles. 
For the exponential three-body force the numerical values of positions and widths 
for the two lowest poles are smaller by about 0.2 MeV and 0.04 MeV, respectively. 
The widths are systematically smaller for the attractive potentials for + and 1~ 
while the 2 + -poles appear to depend on the individual case. 

The widths of these S-matrix poles depend rather sensitively on their energies, 
which are of the same order as the height of the corresponding effective radial bar- 
riers, see Fig. 4. For these states with energies about 1 MeV, any width above 0.4 
MeV corresponds to a smooth structure in the cross sections. Thus, even though 
the three-body interactions only amount to a fine tuning of the energies, the conse- 
quences for the presence and subsequent observation of continuum structures might 
be substantial. 

The low-lying S'-matrix poles seem to be rather close-lying. By a sufficiently 
large additional artificial three-body attraction they move down towards threshold 
and become eventually bound states. Their apparent energies and widths depend 
rather sensitively on the boundary condition introduced when the wave functions 
are matched to the Hankel functions at a given (large) distance p m ax- In general 
the poles move towards the origin until converged with increasing p max . Especially 
the widths are often sensitive to the matching point. They systematically decrease 
or sometimes remain constant with increasing p ma x- 

For the 2 + -resonance the sensitivity is very small when p max is larger than 40 
fm. Most of the other poles, however, require larger p max indicating that they are 
somewhat more related to the larger distances in hyperradius. This might be a 
reminiscence of the Efimov effect, where the bound states are pushed up into the 
continuum, but still with a relatively low-lying and dense energy spectrum. The Efi- 
mov effect would arise as the consequence of very low-lying two-body virtual s-states 
in the neutron-core and the neutron-neutron subsystems. The actual parameters 
give energies of about -200 keV for these virtual s-states^]. 

The main part of the radial wave function is determined by the angular mo- 
mentum dependent effective potential corresponding to the lowest A. We show 
these potentials for ± , 1 ± and 2 ± in Fig. 4, where we also exhibit the parts 
remaining after removal of the generalized centrifugal barrier terms, i.e. h 2 (K + 
3/2)(K + 5/2)/(2mp 2 ), where K(K + 4) = X(p = oo) is the corresponding asymp- 
totic hyperspherical eigenvalue. As for two-body systems this remaining part is 
more revealing than the total potential, which could be repulsive for all p and still 
produce a resonance provided a sufficiently strong attractive pocket is present in 
this "non-centrifugal" part of the potential. 

The pocket in the effective radial potential is absent for angular momentum 
Q — , 1^ and 2~ . The pocket for 2 + , which definitely produces a narrow resonance 
at about 1 MeV, is slightly less pronounced than for + , where a bound state at 
about 1 MeV is present. All the lowest effective potentials are attractive without 
the "centrifugal barrier" and therefore they could give rise to resonances. 

The observables are related to real values of the energy. We therefore solve 
the radial equations on the real axis for energies corresponding both to the real 
values of the S'-matrix pole and to values away from this pole. We show in Fig. 
5 the absolute values of these wave functions for l - and 2 + . The peak at small 
distance is narrower and much more pronounced for 2 + than for 1~. Still for 1~, 
a substantial amount of strength is present between 5 fm and 20 fm whenever the 
energy is within the width of the S-matrix pole. Outside the widths of all the poles 
the wave functions appear with very little probability at distances below 15 fm, see 
the curve for 1~ with the energy 1.4 MeV. 

We have 2 + as a pronounced resonance and 1 _ which only shows up as a much 

1 We shall use negative values corresponding to the energies of the S-matrix poles. 
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Figure 5: The absolute values of the radial wave functions for 6 He (n+n+a) after 
diagonalization as functions of p for real energies in intervals around the real parts 
of the S-matrix poles E r = 0.95 MeV for 1" and E r = 0.845 MeV for 2+. Only the 
dominating component in the full computation is shown. The interactions are the 
same as in Fig. 4. 



smaller and broader peak in Fig. 5. These wave functions reflect the effects on 
the real axis of the properties of the corresponding complex S-matrix poles. For 
other energies and angular momenta a similar picture is found. The small distance 
enhancements are obtained whenever the energy is within the width from the energy 
of an S'-matrix pole. These poles therefore produce observable effects. However, 
the size of the effects depends on both the detailed properties of the poles and the 
precise definition of the observable. It is also clear from Fig. 5 that contributions 
from S'-matrix poles may continuosly vary from substantial to vanishing small. 

3.3 Strength functions and Coulomb cross sections 

The continuum structure can be investigated by electromagnetic and/or nuclear ex- 
citations from the ground state. These transitions are described by observables such 
as the multipole strength functions. The lowest and most important three of these 
are shown in Fig. 6 as functions of energy both for plane waves and for the proper 
continuum wave functions. The curves are normalized by their respective sum rule 
values and each of them would therefore after integration over all energies give 1. 
The strengths below 10 MeV are 97%, 78% and 83% for A = 0, 1,2, respectively 
and for the corresponding plane waves we obtain the somewhat smaller values 94%, 
60% and 81% obtained. 

The influence of the final state interaction is directly reflected in deviations 
from the broader plane wave distributions. In general we always must have a rise 
from zero to a maximum and a fall off towards zero at large energy. Especially 
pronounced peak structure as observed for A = 2 is the signature of a resonance, 
which in this case reflects the well known 2 + state at 0.82 ± 0.025 MeV of width 
0.113 ± 0.020 MeV (49|, which in this computation appears at 0.82 MeV with the 
width 0.093 MeV. 

For 1~ a peak and a shoulder appears at about 0.95 MeV and 1.8 MeV. This 
1 _ enhancement at low energy arises from the two overlapping S-matrix poles seen 
in Table 1, see also J2(|. The enhancement almost coincides in energy with the 
dominating 2 + -peak and consequently it must be harder to detect experimentally. 
The nuclear + strength function resembles the plane wave result more than the 
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Figure 6: The strength functions, dB EX /dE = £J(nJ 7r ||M(£;A)||0 + )| 2 for 
6 He (n+n+a) as functions of energy for transitions from the ground state to + 
(dotted), 1~ (solid) and 2 + (dashed) excited continuum states. The operator is 
M(EX,fi) = p 2 and 5D i=1 eZir^Yx^iji) respectively for A = and 1,2. The units 
are the corresponding sum rule values (0 + |p 4 |0 + ) — (0 + |/9 2 |0 + ) 2 for A = and 
e 2 Z 2 (2A + l)(0+|r 2A |0+)/(47r) for A = 1,2, where eZ a is the 4 He-charge and r a is 
the He-distance from the 6 He center of mass. The interactions are the same as in 
Fig. 4. The smooth curves (smaller at small distance) correspond to plane waves 
for the continuum states. 

higher multipoles reflecting broader underlying structures where the poles have 
larger widths if present at all. 

The differential Coulomb dissociation cross section is now computed by multi- 
plication of strength functions and virtual photon spectra. The results are shown 
in Fig. 7 for Pb and Cu targets for both dipole and quadrupole excitations. As 
expected the dipole contribution has a width of about 2 MeV and it is by far dom- 
inating in absolute size. The quadrupole distribution is much smaller, but strongly 
peaked at the resonance energy. The target dependence vary with the square of the 
target charge as seen from eq. (p0|) . Both potentials as well as different three-body 
interactions give similar, but distinguishable results as seen in the right hand side 
of Fig. 7. The major differences arise from the difference in the ground state struc- 
ture, in particular the larger spatial extension found for the exponential three-body 
force. With the same ground state wave function almost identical strength functions 
would appear. 

Previous computations of l~-strcngth functions reported peaks at about 2.5 
MeV and shoulders at about 6 MeV @ ||, [11 & Tne P re sent 1 "-strength 
function differs substantially with much more strength at low energies indicating 
contributions from larger distances. Unfortunately corresponding experiments are 
so far not available for 6 He. 

The low-energy enhancement of the dipole strength function move strength to- 
wards energies with larger values of the number of virtual photons. The total 
Coulomb dissociation cross section is therefore larger than that obtained with plane 
waves in the final state. It is also necessarily large compared to analogous cross sec- 
tions for ordinary nuclei, again due to a relatively large low-energy enhancement. 
This is explained physically as the result of the weakly bound neutrons easily sep- 
arated by a small Coulomb disturbance. 
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Figure 7: The dipole and quadrupole (inset) contribution to the differential 
Coulomb dissociation cross section of 6 He (n+n+a) at 800 MeV/nucleon as function 
of the three-body energy for a Cu- and a Pb-target. The Cu-results are multiplied 
by (82/29) 2 to remove the dominating, but trivial overall charge scaling. The en- 
ergy of the 6 He-beam is 800 MeV/nucleon. The interaction is the same as in Fig. 4. 
On the right hand side is shown the strength functions dBEi/dE for the attractive 
potential in eq.(|55|) (dashed) and for the repulsive s-wave potential in eq.(Q) (solid 
and dot-dashed) with gaussian (G) and exponential (E) three-body interactions. 
The parameters are given in section 3. 



The total Coulomb dissociation cross section is simply obtained by integrating 
the differential cross section over energy. The quadrupole contribution amounts here 
to about 0.5% and the total cross section is 373 mb and 54 mb for the two targets 
and the beam energy of 800 MeV/A. This is in agreement with the experimental 
extrapolation of J5l| and the calculated values in |22| while somewhat larger than 
computed in piq ] . This rather favorable comparison supports the three-body model 
with a substantial 1~ low-energy enhancement. However, the enhancement is not 
in itself proof of the presence of a low-lying three-body dipole resonance. Any 
attraction would produce more strength at low energies. A resonance needs more 
than marginal attraction. Furthermore, the enhancement could be due to relatively 
strong underlying two-body structures. In the present case the S'-matrix pole at 
about 1 MeV could indicate a resonance, which however overlaps the next pole at 
about 1.5 MeV. This in turn results in the relatively weak peak in the wave function 
at small distances in Figs. 5 and 6. 

The total Coulomb dissociation cross section is shown in Fig. 8 as function of 
beam energy. Both dipole and quadrupole contributions are shown. The experi- 
mental point from |51 has very large error bars and therefore not surprisingly in 
agreement with the computations. At high energy we find a slightly decreasing 
function with values about a few hundred mb. At energies below about 10 MeV/A 
the cross section drops dramatically with decreasing energy. Although the approx- 
imations are dubious at these energies the precise behavior of this rapid change of 
the cross section should be sensitive to the details of the halo structure. 



4 The 11 Li-system as n+n+ 9 Li 

The a-particle has spin zero and the neutron-a system has a low-lying p-resonance. 
Consequently 6 He has spin zero and consists mainly of ^-configurations of the 
neutron-a relative wave function. For n Li the spin and parity is | as for the 9 Li- 
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Figure 8: The total Coulomb dissociation cross section as function of the laboratory 
energy for a 6 He-beam colliding with a 208 Pb-target. Both the quadrupole part and 
the total cross sections are shown. The plane wave results are shown as dashed 
curves. The interaction is the same as in Fig. 4. The experimental point is from 

& 



core of the three-body system. Furthermore the neutron- 9 Li system has apparently 
a virtual s-states at about -200 keV and somewhat higher-lying p-resonances at 



about 600 keV |13, ^3|. Consequently n Li is expected to consist of dominating 
s 2 -configurations of the relative neutron- 9 Li wave function. The details about the 
two-body subsystems are not accurately known, but already the possibility of low- 
lying s-states is interesting, since the conditions for the Efimov effect then nearly 
are fulfilled. In this section we shall therefore use the knowledge obtained from the 
simpler and better known 6 He and predict the more uncertain properties of Li. 



4.1 Interactions and numerical details 

The neutron-neutron interaction is the same as used in the previous computations 
for 6 He. The neutron-core effective interaction often assumes zero spin for both 9 Li 
and n Li although the correct spins are | for both nuclei. The spin-orbit term I ■ s„ 
for a neutron in the relative motion around a spin-zero core is used although the 
natural generalization, which in fact has been used previously || would be of 
the form i ■ (s„ + s c ) for finite core spin s c ^ 0. 

However, in our case of nuclear clusters the Pauli principle must be treated in 
one way or another. This can be achieved either by a large repulsion in forbidden 
states or by omitting or projecting out the forbidden configurations. These different 
approximations assume that the forbidden states can be identified and preferentially 
expressed in terms of the three-body wave functions obtained as solutions for the 
corresponding neutron-core potential. Thus an effective two-body potential is much 
easier to apply in three-body computations when its symmetries, quantum numbers 
and the related eigenfunctions are expressed by the quantities used in calculations 
of the nucleonic motion inside the core. We shall therefore use the mean-field spin- 
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orbit form I • s n , where the spin of the core does not enter. 

In the present work we shall deal with the Pauli principle in several ways and 
then compare the results. The neutrons in the core occupy the lowest Sj/2 and P3/2- 
states. The occupied p 3 / 2 -state is avoided in the three-body computation by using 
a sufficiently large repulsive two-body potential. The neutron- 9 Li s-wave potential 
is either shallow without any bound states or deeper with one bound state but 
the same scattering length and effective range. In the latter case the forbidden 
three-body configuration is excluded in the calculations || . 

Another qualitative difference from the zero core-spin computations is the two 
possible couplings of the spins of the neutron and the core. For the neutron- 9 Li 
the total spin can then be 1 or 2. In general we therefore also include a spin-spin 
potential term to differentiate between these two spin-couplings for each orbital 
angular momentum state. Such spin-splitting terms are most likely present due to 
the strong spin dependence of the underlying basic interaction and consequently 
hard to ignore. 

For finite core spin the interactions corresponding to a shallow s-wave potential 

are 

y^c =Q) = (-7-28 - 0.31s„ • s c ) exp(-r 2 /2.55 2 ) 
v nc =1) = ( 18 -25 + 1.47s„ • s c + • s„) cxp(-r 2 /2.55 2 ) . (56) 

The two s-wave scattering lengths and effective ranges are (7.65 fin, 4.53 fm) and 
(10.88 fm, 4.77 fm) corresponding to virtual s-states at -0.247 MeV and -0.140 
MeV for the total spin of 1 and 2, respectively. The energies and widths of the P1/2- 
resonances defined as a poles of the S'-matrix are E(p 1 / 2 ) — 0.75 MeV, T(p 1 / 2 ) — 
0.87 MeV and E(p 1/2 ) = 1.60 MeV and T(p 1/2 ) = 3.74 MeV, respectively for spin 
1 and 2. In all cases the high- lying p 3 / 2 -resonance is not contributing in the three- 
body calculations and the Pauli blocking by the core neutrons are simulated in this 
way. 

In addition to the two-body potentials a diagonal three-body force could be 
introduced for fine tuning as for 6 He. However, the idea of using the three-body force 
is to include effects beyond those accounted for by the two-body interactions and too 
imprecise and too little information is available about this neutron-core system. It 
is therefore at the moment as reasonable to adjust the two-body interaction instead 
of adding another uncertainty at this level. 

The choice of interaction parameters is dictated by the knowledge of 11 Li and 
the accumulating information about the structure of 10 Li, i.e. a p-resonance at 
about 0.6 MeV, a low-lying virtual s-state and a small spin splitting of these states 
|l3| , [53|]. We obtain a three-body energy of about -300 keV reproducing the 11 Li- 
binding energy of 295 ± 35 keV with the corresponding root mean square radius 
of 3.34 fm. Furthermore, the calculated fragment momentum distributions in n Li 
break-up reactions also compare rather well with measured values Q. Then the n Li 
ground state wave function has about 80% and 20% of s 2 and ^-configurations, 
respectively. 

We shall use this "realistic" interaction in the investigation of the continuum 
properties of 11 Li. All possible s- and p- waves are included. When the large-distance 
asymptotic behavior is reached the solutions are obtained from eqs.(p4[) and (p7|). 
The radial equations are integrated from zero to about 200 fm. Further arguments 
for these numerical choices can be found in [ p6[ . 

4.2 Solutions and S'-matrix poles 

The angular eigenvalues for | , | and | are shown in Fig. 9 together with the 
asymptotic behavior obtained from the analytical expressions. For the ground state 
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Figure 9: The lowest angular eigenvalues A n for n Li (n+n+ 9 Li) as functions of p for 
angular momentum = 5 , 1 , f^. The dot-dashed lines are the large distance 
asymptotic behavior from eq. ( p(i| ) . The neutron-neutron interaction is from || and 
the neutron- 9 Li interaction is given in eq.(^6|). Maximum hyperspherical quantum 
numbers up to about 100 are used. 
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Table 2: The real and imaginary values (E r ,T) (in MeV) of the lowest S-matrix 
poles E — E r — iT/2 for 11 Li for various spins and parities J w . The excitation 
energy E* = E r + 0.305 MeV. The interactions for the upper part of the table are 
the same as in Fig. 1. The middle of the table contains the results for a model with 
s c = and the same average positions of the lowest neutron-core resonances, i.e. 
an average energy of the Si/ 2 virtual state and the p^-resonance at -0.18 MeV and 
1.22 MeV, respectively. The lower part of the table is for a with one Pauli forbidden 
s-state. The si/ 2 virtual state is at -0.18 MeV and the pi/ 2 neutron-core resonance 
is at 0.50 MeV. 





E r r 


E r r 


E r r 


E r r 


E r r 


1- 

2 

3 _ 
2 

5~ 
2 


-0.305 


0.89 0.43 


1.37 0.51 
1.41 0.56 
1.36 0.49 


1.56 0.56 
1.60 0.61 
1.60 0.68 


1.98 0.65 
2.03 0.68 
2.01 0.72 


1 + 
2 

3 + 

1 + 
2 


0.65 0.35 
0.68 0.33 
0.68 0.37 


0.88 0.33 


1.28 0.48 
1.33 0.50 
1.36 0.55 


1.74 0.64 
1.77 0.63 
1.74 0.64 


1.95 0.68 
2.08 0.71 
2.11 0.84 


0+ 

1+ 


-0.307 


1.00 0.37 


1.35 0.45 
1.40 0.59 


1.62 0.61 
1.59 0.63 


1.96 0.92 
2.02 0.81 


0" 

1- 


0.64 0.31 


0.92 0.39 


1.25 0.51 
1.46 0.53 


1.82 0.62 
1.76 0.59 


2.02 0.65 
2.08 0.67 


0+ 

1+ 


-0.306 


1.00 0.31 


1.40 0.41 
1.37 0.49 


1.64 0.56 
1.62 0.58 


1.99 0.88 
2.03 0.75 


o- 
1- 


0.64 0.27 


0.88 0.30 


1.29 0.43 
1.50 0.39 


1.76 0.53 
1.82 0.60 


1.96 0.69 
2.00 0.71 



with | we find the two lowest levels very similar to the spectrum for zero core-spin. 
However, now a series of additional higher-lying levels appear. They arise from the 
broken symmetries due to the finite core spin I. 

The spectra contain almost identical levels for | + and I + as well as for ^ , 

| and | . A number of additional levels are furthermore present for | . The 
corresponding degeneracy is due to the weak neutron-core spin-splitting potential. 
It can be explained by coupling two neutrons in Si> 2 and Pi/ 2 neutron-core states 

to 0^ and 1 — which in turn, coupled to the | from the 9 Li-core, results in sets of 

nearly degenerate ^ and | + , | + -states. The lowest | , | -states arise from 
couplings of higher orbitals. 

For the lowest spins and parities we give in Table 2 the lowest resonance energies 
and the related widths obtained as S-matrix poles by the complex energy method. 
We show the results for core-spins zero with shallow and deep potentials and for the 
spin | . The spin-zero core approximation show a low-lying and relatively narrow 
S-matrix pole for 1~, 0~ and perhaps also for + while we find nothing similar for 
1 + . Higher-lying and broader poles are found for all angular momenta and parities. 
The same structure is found for the deep potential. The shallow and deep potentials 
also quantitatively give very similar results. 

With the correct finite core spin the symmetries are broken. We recognize the 
three times nearly degenerate l~-pole at about 0.65 MeV with a width of about 
0.35 MeV. We also find degenerate 0^ S'-matrix poles at 0.89 MeV with widths of 
0.33 MeV and 0.43 MeV. More discussion about these poles can be found in p(|. 

The relatively large number of S'-matrix poles could be due to the Efimov effect, 
which occurs when the scattering lengths are much larger than the range of the 
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interactions |52], With increasing scattering lengths, the infinitely many poles 
of the three-body S-matrix move towards the point £'=0. For very large but finite 
scattering lengths a number of poles must already appear close to zero. These poles 
originate from the long distance tail of the effective potential (oc X^i=i 3 a iPjkP~ 3 > 
where Oj is the scattering length of the i-th subsystem) and they are not sensitive to 
the details of the interactions. Since there are no confining barriers for these poles, 
their corresponding widths must be rather large. 

For 11 Li the Efimov condition is almost fulfilled, since a nn fj,~^ +2a cn /j,~ r j; «50 fm. 
This must necessarily result in a number of broad S-matrix poles near the E=0 
point. For 6 He there is a low- lying two-body p3/2-resonance while for 11 Li there 
is instead an s^-virtual state. The latter case is therefore closer to the Efimov 
conditions and a larger number of low-lying three-body S-matrix poles could be 
expected. 

An indication of the properties and origin of these S'-matrix poles is obtained 
by matching the radial wave functions at various decreasing values of p m ax- The 
poles move towards larger absolute values of the energies for decreasing matching 
radii. The imaginary values stay almost constant for the lowest + and l~-states 
indicating resonance-like structures. Both imaginary and real values increase for 
the other poles. 

The lowest effective adiabatic potentials determine the radial wave function and 
the energy of possible bound states. They are shown in Fig. 10 for various spins 
and parities. In all cases we find attractive potentials around 50 MeV deep after 
removal of the repulsive centrifugal barriers. Resonances are therefore possible in 
all these channels. With the centrifugal barriers all the potentials are still attractive 
except those corresponding to h and | . The ground state of = | exhibits 
the largest attractive pocket and no barrier for the lowest adiabatic potential and 
a barrier height of 1.7 MeV at about 7 fm for the second potential. The 1~ excited 
states all have attractive pockets as well as repulsive barriers of about 0.6-0.9 MeV 
for p between 10 and 15 fm. 

Compared to + of 6 He we have now a less attractive but broader | -potential 
corresponding to the ground state quantum numbers, see Fig. 4. The 1~ ex- 
cited states for 6 He have no attractive pocket while it is substantial for n Li in 
agreement with the calculated low-lying S'-matrix poles, which appear around the 
barrier height. 

The wave functions corresponding to real energies around the real part of the 
pole energy 0.68 MeV are shown in Fig. 11. For energies below 0.68 MeV and 
outside its width the peak moves to larger distances, but remains of comparable size. 
For the energy 0.8 MeV we find a similar peak at a slightly smaller distance. This 
peak can be viewed as the combined effect on the real axis of the two overlapping 
pole structures at 0.68 MeV and 0.88 MeV, see Table 2. For 1.0 MeV, respectively 
within and outside the widths of the poles at 0.88 MeV and 1.33 MeV, the peak has 
decreased and moved to a smaller distance. None of all these peaks are pronounced 
in comparison with the next peaks of the same wave function. Thus strong l - - 
resonance structures are not obtained. 

4.3 Strength functions and Coulomb cross sections 

The dominating dipole term in electromagnetic excitations can excite the ground 
state to continuum states of J n = ^ + , § + \ § + while the nuclear monopole excitation 

only produce | -states. The corresponding calculated strength functions are shown 
in Fig. 12 together with the results obtained by using plane waves for the continuum 
wave functions. The monopole strength resembles the results of the plane wave 
computation in agreement with the lack of low-lying S-matrix poles below 0.8 MeV, 
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Figure 10: The total effective diagonal radial potential (solid curves) as functions of 
p corresponding to the four lowest A's for J v = | , | , | for n Li (n+n+ 9 Li). The 
dashed curves are the part remaining after removal of the generalized centrifugal 
terms, i.e. h 2 (K + 3/2) (A" + 5/2)/(2mp 2 ), where K(K + 4) = A(p = oo) is the 
corresponding asymptotic hyperspherical eigenvalue. The interactions are as in 
Fig. 9. The insets show the details of the lowest potentials. 
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Figure 11: The absolute values of the radial wave functions for 11 Li (n+n+ 9 Li) as 
functions of p for real energies in an interval around the real part of the S-matrix 

o -f- 

polc E r = 0.68 MeV for | . The interactions are the same as in Fig. 9. 




Figure 12: The strength functions dB EX /dE oc | (nJ*\ \M(E\)\ \ J* °) | 2 , 
M(E\,n) = X)i=i e Zi r iY\n(ri), for xl Li (n+n+ 9 Li) as function of energy for tran- 
sitions from the ground state via + to | (left hand side), via 1 _ (right hand 

side) to | + and | + excited continuum states. The smooth curves (smaller at 
small distance) correspond to plane waves for the continuum states. The curves are 
normalized to the corresponding sum rule values given in Fig. 6. The interactions 
are the same as in Fig. 9. 
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Figure 13: The strength functions dB EX /dE = J2 n j* I (nJ w \\M(EX)\\ Jq )| 2 , 

M(EX,n) = Ei=i eZirfYx^ih), for A = 1 for 11 Li (n+n+ 9 Li) as function of en- 
ergy. The smooth curves correspond to plane waves for the continuum states. The 
interactions are the same as in Table 2 for s c = and Fig. 9 for s c = | except 
for the J JH curve obtained with the potential in [ |45| . The experimental points are 
from Jl3| , |l4| , |l5| ]. The arbitrary units in Jl5| arc normalized to our sum rule value, 
while the absolute data from |l^, [l4| arc left unchanged. 



see Table 2. The dipole strengths are almost proportional to the statistical weights 
of (2 J + 1) and all of them are enhanced significantly above the plane wave results 
at low energies. This enhancement overlaps with the position of the lowest l~-poles 
in Table 2. 

The total dipole strength function, where the contributions from all J rir in the 
continuum are added, is in Fig. 13 compared to the zero core approximation and 
the three available measured distributions. The plane wave result are the same for 
zero and finite core spin, because the ground state essentially is unchanged. The 
interaction with zero core-spin gives a distribution shifted about 100 keV towards 
lower energy compared to the result for the realistic full computation. A lower and 
broader peak is obtained for the potential from 45 where the p 2 -content of the 
three-body wave function is very small. The low-lying l~-poles around 0.65 MeV 
enhance the strength functions at low energies compared to the plane wave results. 

The computed strength functions substantially exceed most of the data points 
| I3| , 0, [l5| in the peak region around 0.55 MeV. (Note however that the data in 
lj| ]l4| contain much less total strength.) A reduction could be achieved with 
higher energy and larger width of the S'-matrix pole, but this would probably only 
be provided by a potential with much too small p 2 -content in the three-body wave 
function. It is also a curious fact that other models provide 1~ strength functions 
substantially closer although not in complete agreement with the data |gl], |55|, |5(| . 
This is of course related to the lack of low-energy 1~ -resonances or S'-matrix poles 
in these computations. 

In this context it is worth pointing out that it is difficult to find the most 
appropriate comparison with the different experimental results in the figure. The 
normalization must be properly chosen and the theoretical results must be folded 
with the distributions (unknown to us) related to the equipment used in the different 
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E lab /nucleon (MeV) 

Figure 14: The dipole approximation (solid curve) to the Coulomb dissociation 
cross section of 11 Li as function of laboratory energy per nucleon for a Pb-target. 
The plane wave results are the dashed curve. The interactions are the same as in 
Fig. 9. The data point is from [ p^[ . 



experiments. 

The strength functions rise from zero to a maximum and then fall off towards 
zero at large energy whereas the virtual photon spectra decrease monotonically with 
excitation energy The low-energy enhancement necessarily implies a larger 

Coulomb cross section, since the dipole is the dominating dissociation mode. We 
show in Fig. 14 computed total cross sections as function of beam energy for a n Li 
projectile on a lead target. The only available data point with relatively large error 



bars is in agreement with the computation 13 , but substantially larger than the 



plane wave result. For 280 MeV/A we find 1458 mb and 1501 mb, respectively for 
the potentials with zero core-spin and with a small p 2 -content of the wave function. 

The cross section decreases at large beam energy due to the similar decrease 
of the virtual photon spectrum. At small energy, where the approximations are 
invalid, we find an increasing cross section due to the low-energy cut-off of the 
strength function by the virtual photon spectrum. At the maximum around 5 MeV 
per nucleon this and other reactions are expected to be sensitive to the details of 
the neutron halo structure. 



5 Summary and conclusions 

Borromean halo systems are almost by definition weakly bound and excited states 
are usually entirely absent. The continuum is therefore easy to excite and unavoid- 
able in descriptions of essentially all reactions involving such particles. The spatial 
extention of the bound state and the small binding energy require rather accurate 
treatment of the large distances. A suitable method was recently developed and 
applied to compute the ground state of halo systems. This method is in the present 
paper extended to apply for the low-energy part of the three-body continuum. 

The Faddeev equations in coordinate space are solved in two steps. First the 
discrete spectrum of the angular part is computed and used as a complete basis set. 
Then the coupled set of radial equations is solved with the appropriate continuum 
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wave boundary conditions. The angular equations for large distances outside the 
short ranges of the potentials are especially simple. However, they are essential and 
therefore treated carefully. 

The three Faddeev components are very useful in the detailed description of the 
particle correlations. We show that only the three s-waves (one in each component) 
for a given total orbital angular momentum couple at large distances. All other par- 
tial waves are decoupled. We therefore solve these much smaller and simpler sets of 
coupled and uncoupled angular equations at large distances. Some of these asymp- 
totic solutions are obtained analytically and others as solutions to trancendcntal 
analytical equations. As a necessary intermediate result we derive a convenient 
expression for the transformation of angular functions between two different Jacobi 
coordinate systems. 

Systems with two identical neutrons and a core of finite spin are specifically 
treated. The continuum spectra of the two Borromean halo nuclei 6 He (n+n+a) 
and n Li (n+n+ 9 Li) are investigated numerically in some detail. Two-body inter- 
actions with and without bound states, but reproducing the observed low-energy 
scattering data are used for 6 He. In addition three-body interactions with several 
radial shapes are added to obtain the measured binding energy. For 11 Li no three- 
body interaction is used, since the two-body interaction is unknown and therefore 
directly parametrized to reproduce anticipated two-body resonances and virtual 
states in addition to the momentum distributions in fragmentation reactions. 

The antisymmetry between the neutrons in the halo and in the core is accounted 
for in two ways. First by using a repulsive or a shallow neutron-core potential 
without bound states. Second by omitting the lowest adiabatic potential arising 
from a more attractive neutron-core potential with one bound state from the set of 
radial equations. The results are compared. 

The adiabatic potentials are decisive for the radial solutions. The lowest po- 
tential in each channel is attractive when the corresponding generalized centrifugal 
barrier is removed. All channels arc therefore potentially able to support resonance- 
like structures. The pocket in the adiabatic radial potential is absent for 6 He for 
angular momentum 0~, and 2 _ and well developed for both + and 2 + . For 11 Li 
the adiabatic potentials for 1~ -excitations all have attractive pockets and effective 
barriers of 0.75 ± 0.15 MeV. Also the + -channel has a well developed attractive 
pocket, but no barrier, for the lowest potential supporting the ground state. The 
second potential for + -excitations is also attractive with a barrier of about 1.7 
MeV. 

We calculate the ^-matrix poles by the complex energy method. The lowest 
of these poles appear slightly above the barriers and their widths are consequently 
relatively large and rather sensitive to fine tuning of the interactions. One exception 
is the known narrow 2 + -resonance in 6 He which is reproduced in the calculation. 
The narrowest low lying poles for 6 He appear for 1~ and 2 _ at about 1 MeV with 
widths of 0.3-0.4 MeV. For n Li they appear for 1~ and ± respectively at about 0.65 
MeV and 0.9 MeV with widths of about 0.35 MeV. The unusually many low-lying 
S'-matrix poles could indicate that the Efimov limit is fairly close. 

We computed the electric excitations from ground to continuum states. The 
strength functions are rather strongly enhanced at low energies due to the low-lying 
S'-matrix poles. The functions extracted from measurements for n Li are apparently 
significantly smaller than our computations. On the other hand the same exper- 
imental information agrees with the computed Coulomb cross section. A proper 
consistent comparison is still lacking. Also for 6 He we obtain enhanced dipole 
strength functions at low energies. Here the experimental information is not avail- 
able, but the observed 1 "-resonance is reproduced almost within the experimental 
uncertainties. We have not attempted to reproduce this somewhat controversial 
resonance more precisely. 
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In conclusion, we have developed a method to solve the three-body problem 
for short-range potentials. The method treats with special care the large distances 
which are essential for the spatially extended halo systems. We investigate the 
continuum spectra for the two halo nuclei 6 He and n Li and find a number of low- 
lying S'-matrix poles. Strength functions are computed and compared with other 
calculations and available experimental data. Various disagreements are pointed 
out and several controversial features are exhibited. 
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Appendix A 

Rotations between different sets of Jacobi coordinates 

We want to "rotate" the wave function from one set of Jacobi coordinates to an- 
other set as defined in eq.jl^). Only the leading order in an expansion in 1/p is 
needed. We must then first express (xj, yj) in terms of (xj, yi). The six-dimensional 
transformation is J5^, |5^| 

Xj = -Xj cos (p k + yi sin <p k , y 3 = -x* sin ip k - y l cos ip k (1) 

where (p k is defined in eq.(p"g|). Defining the angle 7 between Xj and Yi by 

cos7= , (2) 

we have the relation between the hyperangles on and aj related to the two coordi- 
nates 

sin aj — cos ifk sin on + sin ip k cos on + 2 cos Lfk sin (p k sm on cos 014 cos 7 . (3) 

We now expand the following function, related to the function in eq.(p"3|), of aj 
in terms of Legendre polynomials P^(cos7): 

. 7Z : . —, j, = y,A t (ai)P£(cos 7 ) , (4) 

sin ( 2a j) sin x aj cos v aj ^ 

.t'j'L. . 2£+l f 0nlj> y Ls> x s( a :>) . 

A e (oi) = —z— I dcosj " v j, PHCOS7) , (5) 

& J sm[2aj) sm * aj cos » aj 

where aj is the function of a, and 7 defined through eq.(A|^). Then changing the 
integration variable from cos 7 to aj, i.e. 

sin(2a 7 ) 2 

a! cos 7 = — roa,- , (6) 
sin(2ai) sin(2v? fe ) 3 w 

we can rewrite eq.(A^) as 

,j 21+1 pr/2-\ v /2- Vh -ai\ 

A e " (a,) = — — j- — r— ; — j- — 7 / da,- 



sin * a , cos^ a a 



P^(cos7(a i ,a J )) (7) 
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Next we use the following identities 



^(0087) - ^-Yffi^M (8) 



N W j Jy sin< otj co/y ctjYffi (tl Xj , Q Vj ) = ^ ftj^_* 1<a 
xN nlll2 sin^ 1 ^cos fe a l P^ +1 / 2A + 1 / 2 (cos(2a))F^(r! Xi ,^) , (9) 

+ 4 + V+ 1)! 2(2rc + l x + l y + 2) \ 1/2 
V r(n + 4 + 3/2)r(n + ^ + 3/2) / K ' 



B. 



W _ , U L+t 1+ t v+ t V(2*+ l)(2*i + l)(2l 2 + 1) 



x(£iOIO|40)(^ 2 OfO|^0) 



47T 

4 4 * 



K2 ty 



(12) 



where Y[*¥ L is defined in eq.(|12|), N n £ 1 £ 2 are normalzation constants for the Jacobi 



polynomials P n 



^1 + 1/2,^2 + 1/2 



(cos(2a)), {} and () are the 6J-symbols and the Clebsch- 



Gordon coefficients and the coefficients 72. Jf^ e (<fik) are the so-called Raynal- 
Revai coefficients |5§ with if = 4 + ^ = 2n + l x + £ 2 . 
Then by combining eqs.(A^), (A§)-(A|^) we obtain 



b nlj>,Ls>s( a 3) 



£ W^*>%* (<Pk) a t cos* a J P^+ 1 /^+i/2 (cos(2a4)) 



4tt 



(13) 



where the expansion coefficients C are given by 

xKlb^&k)^ a, cos" a l P^ +1 /2/ 2 +i/2 (cos(2ai)) ^i^ , (14) 



Finally we have therefore the desired expression for eq.([13j 



j-> 33 ^ ^ y 



sin(2aj) 

X — r 

sin(2aj) 



dQ Xi dQ yi 



Yi x e y L (^Xi , Qyi) 



Yf, g, L (fl Xj , H Vj ) — C// (aj) 



(15) 



This completes the general derivation of the expression for the transformation of 
the wave function from one set of Jacobi coordinates to another. 
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The large-distance expansion is now found by first approximating eq.(Aph for 
large p and small 014 as ay — <pk + oti cos 7. Then by using on << 1 in eq.(AM) we 
obtain 



sin(2(y9fc) sin^ x tpf. coa e « (fk 



a/' (oi = o) = - - ;.v — — — fo, (16) 



which through eqs.(A14) and (A|l2|) implies that 



Cf} y \*i = 0) = ^0^(0, = 0)^P„ 1 / 2 ' i+1 / 2 (l)^ ^ £ („ fc «£ 
iv ^^ sm(2^fe) sin*- V3fe cos"" ^fc 

1/2 I 1/2 

For 4 = 0,^ = L,n= we can simplify the expression by using P (1) = 

1, ^oiUo-E/^fc) = ( — 1) L cos^'C'Pfc)- We then obtain eq. (|l7|) . The angular wave 
functions corresponding to non-zero i' x rapidly approach zero for large p- values and 
the their contributions are therefore here assumed to be zero. 

Thus for short range interactions for large p and therefore also for small a only 
l Xi = components receive contributions from the rotated wave functions from 
the other Faddeev components. All partial waves with l Xi > can then be solved 
independently. The remaining three components with l Xi = Q,i = 1,2,3 must be 
solved as a set of coupled equations. 
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